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Abstract. There has been remarkable progress recently in both observational and theoretical 
studies of galaxy formation and evolution. Largely due to a combination of deep Hubble Space 
Qs ■ Telescope (HST) imaging, Keck spectroscopy, and COBE far-IR background measurements, new 

^\ I constraints have emerged on the emission history of the galaxy population as a whole. The global 

0^ ■ ultraviolet, optical, near- and far-IR photometric properties of the universe as a function of cosmic 

time cannot be reproduced by a simple stellar evolution model defined by a constant (comoving) 
,^ . star-formation density and a universal (Salpeter) initial mass function, and require instead a 

r^ ' substantial increase in the stellar birthrate with lookback time. While the bulk of the stars 

present today appears to have formed relatively recently, the existence of a decline in the star- 
^sO [ formation density above z fa 2 remains uncertain. The study of the transition from the cosmic 

'dark age' to an ionized universe populated with luminous sources can shed new light on the 
star formation activity at high redshifts, and promises answers to some fundamental questions 
^ . on the formation of cosmic structures. If stellar sources are responsible for photoionizing the 

OO ' intergalactic medium at z « 5, the rate of star formation at this epoch must be comparable or 

OA ■ greater than the one inferred from optical observations of galaxies at z « 3. A population of 

^^ ' quasars at z < 2 could make a significant contribution to the extragalactic background light if 

^«: . dust-obscured accretion onto supermassive black holes is an efficient process. 

i^ ■ INTRODUCTION 



There is little doubt that the last few years have been very exciting times in galaxy 



O 

H , 

c/3 ! formation and evolution studies. The remarkable progress in our understanding of faint 

^ ' galaxy data made possible by the combination of Hubble Space Telescope {HST) deep imaging 
and ground-based spectroscopy has permitted to shed new light on the evolution of the 
^ . stellar birthrate in the universe, to identify the epoch 1 < z <2 where most of the optical 

H \ extragalactic background light was produced, and to set important contraints on galaxy 
evolution scenarios. The explosion in the quantity of information available on the high- 
redshift universe at optical wavelengths has been complemented by the detection of the 
far-IR/sub-mm background by DIRBE and FIRAS onboard the COBE satellite, and by 
theoretical progress made in understanding how cosmic structure forms from initial density 
fluctuations [41]. The IR data have revealed the 'optically-hidden' side of galaxy formation, 
and shown that a significant fraction of the energy released by stellar nucleosynthesis is 
re-emitted as thermal radiation by dust. The underlying goal of all these efforts is to 
understand the growth of cosmic structures, the internal properties of galaxies and their 
evolution, the mechanisms that shaped Bubble's morphological sequence, and ultimately to 
map the transition from the cosmic 'dark age' to a ionized universe populated with luminous 
sources. While one of the important questions recently emerged is the nature (starbursts 



^^ To appear in Physica Scripta, Proceedings of the Nobel Symposium, Particle Physics and the Universe 
(Enkoping, Sweden, August 20-25, 1998). 
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FIGURE 1. Differential galaxy number counts per square degree as a function of apparent magnitude in 
four bandpasses from near-UV to near-IR. The sources of the data points are indicated in each panel (see [44] 
for the complete reference list). Note the decrease of the logarithmic slope d\ogN/dm at faint magnitudes. 

or active galactic nuclei?) and redshift distribution of the ultraluminous sub-mm sources 
discovered by SCUBA, of perhaps equal interest is the possible existence of a large population 
of faint galaxies still undetected at high redshifts, as the color-selected ground-based and 
Hubble Deep Field (HDF) samples include only the brightest and bluest star-forming objects. 
In any hierarchical clustering ('bottom-up') scenario (the cold dark matter model being the 
best studied example), subgalactic structures are the first non-linearities to form. High- 2; 
dwarf galaxies and/or mini-quasars (i.e. an early generation of stars and accreting black 
holes in dark matter halos with circular velocities Vc ~ 50 kms~^) may then be one of the 
main source of UV photons and heavy elements at early epochs. 

In this talk I will focus on some of the open issues and controversies surrounding our 
present understanding of the history of the conversion of cold gas into stars within galaxies, 
and of the evolution of luminous sources in the universe. An Einstein-deSitter (EdS) universe 
[Qm = 1, ^A = 0) with h = i/o/lOO kms~^ Mpc^^ = 0.5 will be adopted in the following. 

COUNTING GALAXIES 



Much observing time has been devoted in the past few years to the problem of the detection 
of galaxies at high redshifts, as it was anticipated that any knowledge of their early luminosity 
and color evolution would set important constraints on the history of structure and star 
formation in the universe. As the best view to date of the optical sky at faint flux levels, the 
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HDF imaging survey has rapidly become a key testing ground for models of galaxy evolution. 
The field, an undistinguished portion of the northen sky at high galactic latitudes (the data 
from a southern deep field are being analyzed as we speak), is essentially a deep core sample 
of the universe, acquired with the HST in a 10-day exposure. With its depth - reaching 
5-(T limiting AB magnitudes of roughly 27.7, 28.6, 29.0, and 28.4 in f/, 5, V, and / ^ - and 
four-filter strategy to provide constraints on the redshift and age distribution of galaxies in 
the image, the HDF has offered the astronomical community the opportunity to study the 
galaxy population in unprecedented detail [57]. 

There are about 3000 galaxies in the HDF, corresponding to 2 x 10® deg~^ down to the 
faint limit of the images. The galaxy counts are shown in Figure 1 in four bandpasses 
centered at roughly 300, 450, 600, and 800 nm. A compilation of existing ground-based 
data is also shown, together with the predictions of no-evolution models, i.e. models in 
which the absolute brightness, volume density, and spectra of galaxies do not change with 
time. In all four bands, the logarithmic slope a of the galaxy number-apparent magnitude 
counts, logA^(m) = am, flattens at faint magnitudes, e.g., from a = 0.45 in the interval 
21 < 5 < 25 to a = 0.17 for 25 < -B < 29. The slope of the galaxy counts is a simple 
cosmological probe of the early history of star formation. The flattening at faint apparent 
magnitudes cannot be due to the reddening of distant sources as their Lyman break gets 
redshifted into the blue passband,^ since the fraction of Lyman-break galaxies at i? ~ 25 is 
only of order 10%. Moreover, an absorption-induced loss of sources could not explain the 
similar flattening of the galaxy counts observed in the V and / bands. Rather, the change 
of slope suggests that the surface density of luminous galaxies declines beyond z ~ 1.5. 

THE BRIGHTNESS OF THE NIGHT SKY 

The extragalactic background light (EBL) is an indicator of the total luminosity of the uni- 
verse. It provides unique information on the evolution of cosmic structures at all epochs, as 
the cumulative emission from galactic systems and active galactic nuclei (AGNs) is expected 
to be recorded in this background. 

The contribution of known galaxies to the optical EBL can be calculated directly by 
integrating the emitted flux times the differential galaxy number counts down to the de- 
tection threshold. The leveling off of the counts is clearly seen in Figure 2, where the 
function i^ = iQ-0-^i™-+^8-^) x N{m) is plotted against apparent magnitude in all bands [44]. 
While counts having a logarithmic slope of a > 0.40 continue to add to the EBL at the 
faintest magnitudes, it appears that the HDF survey has achieved the sensitivity to capture 
the bulk of the extragalactic light from discrete sources (an extrapolation of the observed 
counts to brighter and/or fainter magnitudes would typically increase the sky brightness by 
less than 20%). To AB = 29 mag, the sky brightness from resolved galaxies in the /-band 
is ~ 2 X 10~^°ergscm~^ s~^ Hz^^ sr^^, increasing roughly as A^ from 2000 to 8000 A. The 



^•' To get a feeling of the depth of this survey, note that AB = 29 mag corresponds to the flux at Earth 

from a 100 Watt light bulb at a distance of 10 million kilometers. 

3) Por galaxies with z > 2 (z > 3.5), the H i Lyman edge shifts into the 300 (450) nm HDF bandpass. 

Neutral hydrogen, which is ubiquitous both within galaxies and in intergalactic space, strongly absorbs 

ultraviolet light, creating a spectral discontinuity that can be used to identify young, high- redshift galaxies 

[50]. 
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FIGURE 2. The contribution of known galaxies to the extragalactic background Ught per magnitude bin 
as a function of U {filled circles), B [open circles), V {filled triangles), I {open squares) and K {filled squares) 
magnitudes. For clarity, the B, V, I and K values have been multiplied by a factor of 2, 8, 20, and 40, 
respectively. 



flattening of tlie number counts lias the interesting consequence that the galaxies that pro- 
duce ~ 60% of the blue EBL have B < 24.5. They are then bright enough to be identified 
in spectroscopic surveys, and are indeed known to have median redshift {z) = 0.6 [32]. The 
quite general conclusion is that there is no evidence in the number-magnitude relation down 
to very faint flux levels for a large amount of star formation at high redshift. Note that 
these considerations do not constrain the rate of starbirth at early epochs, only the total 
(integrated over cosmic time) amount of stars - hence background light - being produced, 
and neglect the effect of dust reddening. 

Figure 3 shows the total optical EBL from known galaxies together with the recent COBE 
results. The value derived by integrating the galaxy counts [44] down to very faint magnitude 
levels [because of the flattening at faint magnitudes of the number-magnitude relation most 
of the contribution to the optical EBL comes from relatively bright galaxies] implies a lower 
limit to the EBL intensity in the 0.3-2.2 fun interval of /opt ~ 12 nWm~^sr~^.^ When 
combined with the FIRAS and DIRBE measurements (/fir ~ 16 nWm~^sr~^ in the 125- 
5000 yum range), this gives an observed EBL intensity in excess of 28 nWm~^sr~^. The 
correction factor needed to account for the residual emission in the 2.2 to 125 /im region 
is probably ^ 2 [11]. (We shall see below how a population of dusty AGNs could make 
a signiflcant contribution to the FIR background.) In the rest of this talk I will adopt 
a conservative reference value for the total EBL intensity associated with star formation 
activity over the entire history of the universe of /ebl = 40 /40 nW m"^ sr~^. 



"*) The direct detection of the optical EBL at 3000, 5500, and 8000 A derived from HST data [3] imphes 
values that are about a factor of two higher than the integrated light from galaxy counts. 
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FIGURE 3. Spectrum of the extragalactic background light as derived from a compilation of ground-based 
and space-based galaxy counts in the U, B, V, I, and X-bands {filled dots), together with the FIRAS 125-5000 
/im (solid and dashed lines) and DIRBE 140 and 240 fim {filled squares) detections [15], [25]. The empty 
squares show the DIRBE points after correction for WIM dust emission [29] . 

MODELING GALAXY EVOLUTION 



In the past few years two different approaches have been widely used to interpret faint 
galaxy data [12]. In the simplest version of what I will call the 'traditional' scheme, a 
one-to-one mapping between galaxies at the present epoch and their distant counterparts 
is assumed: one starts from the local measurements of the distribution of galaxies as a 
function of luminosity and Hubble type and models their photometric evolution assuming 
some redshift of formation and a set of parameterized star formation histories [52]. These, 
together with an initial mass function (IMF) and a cosmological model, are then adjusted to 
match the observed number counts, colors, and redshift distributions. Beyond the intrinsic 
simplicity of assuming a well defined collapse epoch and pure-luminosity evolution thereafter, 
the main advantage of this kind of approach is that it can easily be made consistent with 
the classical view that ellipticals and spiral galaxy bulges (both redder than spiral disks 
and containing less gas) formed early in a single burst of duration 1 Gyr or less [2] . Spiral 
galaxies, by contrast, are characterized by a slower metabolism, i.e. star formation histories 
that extend to the present epoch. In these models, typically, much of the action happens at 
high-redshifts. 

A more physically motivated way to interpret the observations is to construct semianalytic 
hierarchical models of galaxy formation and evolution [56]. Here, one starts ab initio from a 
power spectrum of primordial density fluctuations, and follows the formation and hierarchical 
merging of the dark matter halos that provide the early seeds for later galaxy formation. 
Baryonic gas gets accreted onto the halos and is shock-heated. Various prescriptions for 
gas cooling, star formation, feedback, and dynamical friction are adopted, and tuned to 
match the statistical properties of both nearby and distant galaxies. In this scenario, there 
is no period when bulges and ellipticals form rapidly as single units and are very bright: 
rather, small objects form first and merge continually to make larger ones. Galaxy do not 
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FIGURE 4. Left: Mean comoving density of star formation as a function of cosmic time. The data points 
with error bars have been inferred from the UV-continuum luminosity densities of [32] {filled dots), [8] {filled 
squares), [37] {filled pentagons) , [54] {empty dot), and [51] {empty square). The dotted line shows the fiducial 
rate, (p,) — 0.054 M0 yr^^ Mpc^'^, required to generate the observed EBL. Right: dust corrected values 
(^1500 = 1-2 mag). The Ha determinations of [17], [53], and [19] {filled triangles), together with the SCUBA 
lower limit [27] {empty pentagon) have been added for comparison. 

evolve as isolated objects, and the rate of interaction was higher in the past. The bulk of 
the galaxy population is predicted to have been assembled quite recently, and most galaxies 
never experience star formation rates in excess of a few solar masses per year. 



STAR FORMATION HISTORY 



Recently, it has become familiar to follow an alternative method, which focuses on the 
emission properties of the galaxy population as a whole. It traces the cosmic evolution with 
redshift of the galaxy luminosity density and offers the prospect of an empirical determi- 
nation of the global star formation history of the universe and IMF of stars independently 
of the merging histories, complex evolutionary phases, and possibly short-lived star forma- 
tion episodes of individual galaxies. The technique relies on two basic properties of stellar 
populations: a) the UV-continuum emission in all but the oldest galaxies is dominated by 
short-lived massive stars, and is therefore a direct measure, for a given IMF and dust content, 
of the instantaneous star formation rate; and b) the rest-frame near-IR light is dominated by 
near-solar mass evolved stars, the progenitors of which make up the bulk of a galaxy's stellar 
mass, and is more sensitive to the past star-formation history than the blue (and UV) light. 
By modeling the "emission history" of the universe at ultraviolet, optical, and near-infrared 
wavelengths from the present epoch to high redshifts, one should be able to shed light on 
some key questions in galaxy formation and evolution studies: Is there a characteristic epoch 
of star and metal formation in galaxies? What fraction of the luminous baryons observed 
today were already locked into galaxies at early epochs? Are high- 2; galaxies obscured by 
dust? Do spheroids form early and rapidly? Is there a universal IMF? 

The comoving volume-averaged history of star formation follows a relatively simple de- 
pendence on redshift. Its latest version, uncorrected for dust extinction, is plotted in Figure 
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4 {left). The measurements are based upon the rest-frame UV luminosity function (at 1500 
and 2800 A), assumed to be from young stellar populations [34]. The prescription for a 'cor- 
rect' de- reddening of these values has been the subject of an ongoing debate. Dust may play 
a role in obscuring the UV continuum of Canada-France Reshift Survey (CFRS, 0.3 < z < 1) 
and Lyman-break {z ~ 3) galaxies, as their colors are too red to be fitted with an evolving 
stellar population and a Salpeter IMF [37]. Figure 4 (right) depicts an extinction-corrected 
version of the same plot. The best-fit cosmic star formation history (shown by the dashed- 
line) with such a universal correction produces a total EBL of 37 nWm~^ sr~^. About 65% 
of this is radiated in the UV+optical+near-IR between 0.1 and 5 /xm; the total amount of 
starlight that is absorbed by dust and reprocessed in the far-IR is 13 nWm~^ sr^^. Because 
of the uncertainties associated with the incompleteness of the data sets, photometric redshift 
technique, dust reddening, and UV-to-SFR conversion, these numbers are only meant to be 
indicative. On the other hand, this very simple model is not in obvious disagreement with 
any of the observations, and is able, in particular, to provide a reasonable estimate of the 
galaxy optical and near-IR luminosity density. 

THE STELLAR BARYON BUDGET 

With the help of some simple stellar population synthesis tools it is possible at this stage 
to make an estimate of the stellar mass density that produced the integrated light observed 
today. The total bolometric luminosity of a simple stellar population (a single generation 
of coeval, chemically homogeneous stars) having mass M can be well approximated by a 
power-law with time for all ages t > 100 Myr, 



L{t) = 1.3 Lq 



-0.8 



MoVlGyr; 

(cf. [6]), where we have assumed solar metallicity and a Salpeter IMF truncated at 0.1 and 
125 Mq. In a stellar system with arbitrary star formation rate per unit cosmological volume, 
p*, the comoving bolometric emissivity at time t is given by the convolution integral 

Pbol(t) = / L{T)p,{t-T)dT. (2) 

Jo 
The total background light observed at Earth (t = tn) is 

/,e. = fr"f^d, (3) 

4:11 Jo 1 + Z 

where the factor {1 + z) at the denominator is lost to cosmic expansion when converting 
from observed to radiated (comoving) luminosity density. From the above equations it is 
easy to derive in a EdS cosmology 



-2-1/ P* \ / ^H 



1.87 



^-"--^^°°^ "'"-"< Mey.-'Mpc-3 ) (13%; j . (4) 

The observations shown in Figure 3 therefore imply a "fiducial" mean star formation density 
of (p*) = 0.054/40 M0 yr"-*^ Mpc"'^. The total stellar mass density observed today is 
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FIGURE 5. Left: Synthetic evolution of the near-IR luminosity density at rest-frame wavelengths of 1.0 
[long-dashed line) and 2.2 /im {short-dashed line). The model assumes a constant star formation rate of 
p* = 0.054 Mq yr"i Mpc"^ (Salpeter IMF). The dotted (2.2 /xm) and dash-dotted (1.0 /xm) curves show the 
emissivity of a simple stellar population with formation redshift Zon = 5, and total mass equal to the mass 
observed in spheroids today [16]. The data points are taken from [32] {filled dots) and [18] {filled triangle). 
Right: Same but in the i3-band. The data points are taken from [32] {filled dots), [13] {empty triangles), 
and [8] {empty squares). 



P*{tH) = {l-R) p*{t)dt ^ 5 X 10^/40 Mq Mpc" 



(5) 



(corresponding to fi* = 0.007 /40), where R is the mass fraction of a generation of stars 
that is returned to the interstellar medium, R ^ 0.3 for a Salpeter IMF. The optical/FIR 
background therefore requires that about 10% of the nucleosynthetic baryons today [5] are 
in the forms of stars and their remnants. The predicted stellar mass-to-blue light ratio is 
[M/ Lb) ~ 5. These values are quite sensitive to the lower-mass cutoff of the IMF, as very- 
low mass stars can contribute significantly to the mass but not to the integrated light of 
the whole stellar population. A lower cutoff of 0.5 M© instead of the 0.1 M© adopted would 
decrease the mass-to-light ratio (and fi*) by a factor of 1.9 for a Salpeter function. 



TWO EXTREME SCENARIOS 

Based on the agreement between the z k. ^ and z k. A luminosity functions at the bright 
end, it has been recently argued [51] that the decline in the luminosity density of faint HDF 
Lyman-break galaxies observed in the same redshift interval [34] may not be real, but simply 
due to sample variance in the HDF. When extinction corrections are applied, the emissivity 
per unit comoving volume due to star formation may then remain essentially flat for all 
redshift z > 1 (see Fig. 4). While this has obvious implications for hierarchical models of 
structure formation, the epoch of first light, and the reionization of the intergalactic medium 
(IGM), it is also interesting to speculate on the possibility of a constant star-formation 
density at all epochs < 2; < 5. Figure 5 shows the time evolution of the blue and near- 
IR rest-frame luminosity density of a stellar population characterized by a Salpeter IMF, 
solar metallicity, and a (constant) star formation rate of p* = 0.054 MQyr"^ Mpc~^ (needed 



to produce the observed EBL). The predicted evolution appears to be a poor match to the 
observations: it overpredicts the local B and K-hand luminosity densities, and underpredicts 
the 1 /xm emissivity at z ~ 1 . 

At the other extreme, we know from stellar population studies that about half of the 
present-day stars are contained in spheroidal systems, i.e. elliptical galaxies and spiral 
galaxy bulges, and that these stars formed early and rapidly. The expected rest-frame blue 
and near-IR emissivity of a simple stellar population with formation redshift Zon = 5 and 
total mass density equal to the mass in spheroids observed today is shown in Figure 5. In 
this model the near-IR and blue emissivities at z = 4 — 5 are comparable with the values 
observed at z = 1. iJST-NICMOS deep observations may be able to test similar scenarios 
for the formation of elliptical galaxies at early times. 

THE MASS DENSITY IN BLACK HOLES 

Recent dynamical evidence indicates that supermassive black holes reside at the center of 
most nearby galaxies. The available data (about 30 objects) show a strong correlation (but 
with a large scatter) between bulge and black hole mass [38], with Mbh = 0.006 Mbuigc as a 
best-fit. The total mass density in spheroids today is r2buigc = 0.0036lo;ooi7 [16], implying a 
mean mass density of dead quasars 

pi^i^ = 1.5+1-0 X 10^ Mq Mpc-^ (6) 

Since the observed energy density from all quasars is equal to the emitted energy divided 
by the average quasar redshift [49], the total contribution to the EBL from accretion onto 
black holes is 

/bh = ^JT^ ^ 18 nWm-'sT'%,,{l + z)-\ (7) 

4:71 (1 + Z) 

where r^o.i is the efficiency for transforming accreted rest-mass energy into radiation (in units 
of 10%). Quasars at z < 2 could then make a significant contribution to the brightness of 
the night sky if dust-obscured accretion onto supermassive black holes is an efficient process 

[21], [14]. 5 

THE END OF THE 'DARK AGES' 

The epoch of reionization marked the end of the 'dark ages' during which the ever-fading 
primordial background radiation cooled below 3000 K and shifted first into the infrared 
and then into the radio. Darkness persisted until early structures collapsed and cooled, 
forming the first stars and quasars that lit the universe up again [45]. The application of 
the Gunn- Peterson constraint on the amount of smoothly distributed neutral material along 
the line of sight to distant objects requires the hydrogen component of the diffuse IGM to 
have been highly ionized hj z ^ 5 [47], and the helium component hj z ^ 2.5 [10]. From 



^> It might be interesting to note in this context that a population of AGNs with strong intrinsic absorption 
(Type II quasars) is actually invoked in many current models for the X-ray background [33], [7]. 



10 




2 3 4 

redshift 



2 3 4 
redshift 



FIGURE 6. Left: comoving space density of bright QSOs as a function of redshift. The data points with 
error bars are taken from [24] (filled dots), [55] (filled squares), [46] (crosses), and [28] (filled pentagon). The 
empty triangles show the space density of radio- loud quasars [26] . Right: comoving emission rate of hydrogen 
Lyman-continuum photons {solid line) from QSOs, compared with the minimum rate {dashed line) which is 
needed to fully ionize a fast recombining (with gas clumping factor C ~ 30) EdS universe with ilf,^^ = 0.02. 
Models based on photoionization by quasar sources appear to fall short at z « 5. The data point shows the 
estimated contribution from star- forming galaxies at z « 3, assuming that the fraction of Lyman continuum 
photons which escapes the galaxy H i layers into the intergalactic medium is 50% (see [35] for details). 



QSO absorption studies we also know that neutral hydrogen at high-z accounts for only a 
small fraction, ~ 10%, of the nucleosynthetic baryons [30]. 

A substantial population of dwarf galaxies having star formation rates < O.SMoyr^^, 
and a space density in excess of that predicted by extrapolating to faint magnitudes the 
best-fit Schechter function, may be expected to form at early times in hierarchical clustering 
models, and has been recently proposed [39], [35] as a possible candidate for photoionizing 
the IGM at these early epochs. Establishing the character of cosmological ionizing sources 
is an efficient way to constrain competing models for structure formation in the universe, 
and to study the collapse and cooling of small mass objects at early epochs. 

The study of the candidate sources of ionization at z = 5 can be simplified by noting that 
the breakthrough epoch (when all radiation sources can see each other in the hydrogen Lyman- 
continuum) occurs much later in the universe than the overlap epoch (when individual ionized 
zones become simply connected and every point in space is exposed to ionizing radiation). 
This implies that at high redshifts the ionization equilibrium is actually determined by the 
instantaneous UV production rate [35]. The fact that the IGM is rather clumpy and still 
optically thick at overlapping, coupled to recent observations of a rapid decline in the space 
density of radio- loud quasars and of a large population of star- forming galaxies at z > 3, 
has some interesting implications for rival ionization scenarios and for the star formation 
activity in the interval < 3 < 2; < 5. 

The existence of a decline in the space density of bright quasars at redshifts beyond 
~ 3 was first suggested by [40], and has been since then the subject of a long-standing 
debate. In recent years, several optical surveys have consistently provided new evidence for 
a turnover in the QSO counts [24] , [55] , [46] , [28] . The interpretation of the drop-off observed 
in optically selected samples is equivocal, however, because of the possible bias introduced 
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by dust obscuration arising from intervening systems. Radio emission, on the other hand, is 
unaffected by dust, and it has recently been shown [48] that the space density of radio- loud 
quasars also decreases strongly for z > 3. This argues that the turnover is indeed real and 
that dust along the line of sight has a minimal effect on optically-selected QSOs. In this 
case the QSO emission rate of hydrogen ionizing photons per unit comoving volume drops 
by a factor of 3 from z = 2.5 to z = 5, as shown in Figure 6. 

Galaxies with ongoing star-formation are another obvious source of Lyman-continuum 
photons. Since the rest-frame UV continuum at 1500 A (redshifted into the visible band for 
a source at 2; ~ 3) is dominated by the same short-lived, massive stars which are responsible 
for the emission of photons shortward of the Lyman edge, the needed conversion factor, about 
one ionizing photon every 10 photons at 1500 A, is fairly insensitive to the assumed IMF 
and independent of the galaxy history for t 3> 10^ yr. Figure 6 (right) shows the estimated 
Lyman-continuum luminosity density of galaxies at z ^ 3. The data point assumes a value 
of /esc = 0.5 for the unknown fraction of ionizing photons which escapes the galaxy H i layers 
into the intergalactic medium. One should note that, while highly reddened galaxies at high 
redshifts would be missed by the Lyman-break color technique (which isolates sources that 
have blue colors in the optical and a sharp drop in the rest-frame UV), it seems unlikely 
that very dusty objects (with /esc <^ 1) would contribute in any significant manner to the 
ionizing metagalactic flux. 

REIONIZATION 

When an isolated point source of ionizing radiation turns on in a neutral medium, the 
ionized volume initially grows in size at a rate fixed by the emission of UV photons, and an 
ionization front separating the H 11 and H i regions propagates into the neutral gas. Most 
photons travel freely in the ionized bubble, and are absorbed in a transition layer. The 
evolution of an expanding H 11 region is governed by the equation 

'"'■-3HV, *>» ''' 



dt uh t 
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where V/ is the proper volume of the ionized zone, A^ion is the number of ionizing pho- 
tons emitted by the central source per unit time, nH is the mean hydrogen density of the 
expanding IGM, H is the Hubble constant, and tree is the hydrogen mean recombination 
timescale, 

tree= [(l + 2x)nHa,,C]-i = 0.3Gyr(^^) (^) '^30'- (9) 

One should point out that the use of a volume-averaged clumping factor, C = (^hii)/"^?!!!' 
in the recombination timescale is only justified when the size of the H 11 region is large 
compared to the scale of the clumping, so that the effect of many clumps (filaments) within 
the ionized volume can be averaged over. The validity of this approximation can be tested 
by numerical simulations (see Figure 7). Across the I- front the degree of ionization changes 
sharply on a distance of the order of the mean free path of an ionizing photon. When tree 
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FIGURE 7. Simulating the reionization of the universe: propagation of an ionization front in a 128'^ 
cosmological density field. A 'mini-quasar' with A^ = 5 x lO'^^ s^^ was turned on at the densest cell, in a 
virialized halo of total mass 1.3 x 10^^ Mq. The box length is 2.4 comoving Mpc. The solid contours give 
the position of the front at 0.15, 0.25, 0.38, and 0.57Myr after the quasar has switched on at z = 7. The 
underlying greyscale image indicates the initial H i density field. (From [1].) 



is much smaller than the Hubble tinie, the growth of the H ii region is slowed down by 
recombinations in the highly inhomogeneous medium, and its evolution can be decoupled 
from the expansion of the universe. 

In analogy with the individual H ii region case, it can be shown that the hydrogen compo- 
nent in a highly inhomogeneous universe is completely reionized when the number of photons 
emitted above 1 ryd in one recombination time equals the mean number of hydrogen atoms. 
At any given epoch there is a critical value for the emission rate of ionizing photons per unit 
cosmological comoving volume. 



A/ic 



riH(O) 



:io 



51.2 „-l 



s"' Mpc"'') C; 



30 



6 



'Vtnh'^V 



0.02 



(10) 



which is independent of the (unknown) previous emission history of the universe: only rates 
above this value will provide enough UV photons to ionize the IGM by that epoch. One 
can then compare our estimate of A/ion to the inferred contribution from QSOs and star- 
forming galaxies. The uncertainty on this critical rate is difficult to estimate, as it depends 
on the dumpiness of the IGM (scaled in the expression above to the value inferred at ^ = 5 
from numerical simulations [20]) and the nucleosynthesis constrained baryon density. The 
evolution of the critical rate as a function of redshift is plotted in Figure 6 {right). While 
A/ion is comparable to the quasar contribution at z ^ 3, there is some indication of a deficit 
of Lyman-continuum photons at z = 5. For bright, massive galaxies to produce enough UV 
radiation at 2; = 5, their space density would have to be comparable to the one observed 
at z ~ 3, with most ionizing photons being able to escape freely from the regions of star 
formation into the IGM. This scenario may be in confiict with direct observations of local 
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starbursts below the Lyman limit showing that at most a few percent of the stellar ionizing 
radiation produced by these luminous sources actually escapes into the IGM [31].^ 

It is interesting to convert the derived value of A/ion into a "minimum" star formation rate 
per unit (comoving) volume, p*: 

p.{z)=UUz) X lO-^^-Ve^c ^ 0.013/-,^ (^) M^yr^^Mpc-^. (11) 

The star formation density given in the equation above is comparable with the value directly 
"observed" (i.e., uncorrected for dust reddening) at z ~ 3 [37]. The conversion factor 
assumes a Salpeter IMF with solar metallicity, and has been computed using a population 
synthesis code [4]. It can be understood by noting that, for each 1 Mq of stars formed, 8% 
goes into massive stars with M > 20Mq that dominate the Lyman-continuum luminosity 
of a stellar population. At the end of the C-burning phase, roughly half of the initial mass 
is converted into helium and carbon, with a mass fraction released as radiation of 0.007. 
About 25% of the energy radiated away goes into ionizing photons of mean energy 20 eV. 
For each 1 Mq of stars formed every year, we then expect 



0.08 X 0.5 X 0.007 x 0.25 x MqC^ 1 



10^^ phot s~^ (12) 



20 eV 1 yr 

to be emitted shortward of 1 ryd. 

CONCLUSIONS 

Recent studies of the volume-averaged history of stellar birth are pointing to an era of 
intense star formation aX z ^ 1 — 1.5. The optical datasets imply that a fraction close to 
65% of the present-day stars was produced at 2; > 1, and only 25% at z > 2. About half of 
the stars observed today would be more than 9 Gyr old, and only 10% would be younger 
than 5 Gyr.^ There is no 'single epoch of galaxy formation': rather, it appears that galaxy 
formation is a gradual process. Numerous uncertainties remain, however, particularly the 
role played by dust in obscuring star- forming objects. Our first glimpse of the history of 
galaxies to 2; ~ 4 leads to the exciting question of what happened before. Substantial 
sources of ultraviolet photons must have been present at z > 5 to keep the universe ionized, 
perhaps low-luminosity quasars [23] or a first generation of stars in dark matter halos with 
virial temperature T^[j- ~ 10^ — 10^ K [42], [22]. Early star formation provides a possible 
explanation for the widespread existence of heavy elements in the Lya forest [9], while 
reionization by QSOs may produce a detectable signal in the radio extragalactic background 
at meter wavelengths [36]. A detailed exploration of such territories must await projected 
facilities like the Next Generation Space Telescope and the Square Kilometer Radio Telescope. 

Acknowledgments 

Support for this work was provided by NASA through ATP grant NAG5-4236. 



^) At z = 3 Lyman-break galaxies radiate into the IGM more ionizing photons than QSOs if /esc ^ 30%. 
^^ Unlike the measured number densities of objects and rates of star formation, the integrated stellar mass 
density does not depend on the assumed cosmological model. 
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